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Abstract 

Gravitational core-mantle coupling may be the cause of the observed variable 
acceleration of the Earth’s rotation on the 1000 year timescale. The idea is that density 
inhomgeneities which randomly come and go in the liquid outer core gravitationally attract 
density inhomogeneities in the mantle and crust, torquing the mantle and changing its 
rotation state. The corresponding torque by the mantle on the core may also explain the 
westward drift of the magnetic field of 0.2° per year. Gravitational core-mantle coupling 
would stochastically affect the rate of change of the Earth’s obliquity by just a few per cent. 
Its contribution to polar wander would only be about 0.5% the presently observed rate. 
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Tidal friction is slowing down the rotation of the Earth, overwhelming a smaller 
positive acceleration from postglacial rebound. According to historical records such as 
ancient and medieval observations of solar eclipses ( 1 , 2 ), the signal from these two secular 
effects appears to be significantly corrupted by a source (or sources) of an episodic or 
perhaps periodic nature operating on a 1000 year timescale (Figure 1). This nonsecular 
source causes accelerations which are the same order-of-magnitude as the tides, ~6 xlO' 22 s' 
2 . Coupling between the liquid outer core of the Earth and the mantle has long been a 
suspected reason for changes in the length-of-day ( 3 , 4 ). Research heretofore has been 
aimed at viscous, electromagnetic, and topographic coupling at the core-mantle boundary 
( 5 ), mostly on the decadal timescale, and at gravitational coupling not related to the problem 
outlined below ( 6 ). 

The present investigation focusses on the gravitational coupling between the density 
anomalies in the convecting liquid outer core and those in the mantle and crust as a possible 
cause for the observed nonsecular acceleration on the millenial timescale. (Hereafter I use 
"mantle" as a shorthand for "mantle and crust.") The basic idea is as follows (Figure 2). 
There are density inhomogeneities caused by blobs circulating in the outer core like the 
blobs in a lava lamp(7,S); thus the outer core's gravitational field is not featureless. 
Moreover, these blobs will form and dissipate somewhat randomly. Thus there will be a 
time variability to the field. These density inhomogeneities will gravitationally attract the 
density anomalies in the mantle. 

The outer core and mantle may not be in a state of gravitational equilibrium if the 
blobs form and dissipate faster than either the outer core or mantle can react to achieve a 
state of zero torque. In other words, the mantle and outer core may always be in a state of 
disequilibrium due to the stochastic nature of "core weather." The outer core and mantle 
will try to align with each other to make the torque between them zero; but before this can 
happen, the gravitational field of the outer core may change, so that the core and mantle 
never catch up with each other. Hence it may be that the outer core is forever torquing the 
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mantle, changing its rotation state and perhaps producing the observed nonsecular 
acceleration and changing the length-of-day. By action-reaction the outer core will also 
change its rotation state. A quantitive estimate of this effect follows. 

The gravitational potential at a point (r,0,A) exterior to the Earth can be expressed in 
the usual spherical harmonic form 

V(r, 0,X) = X^os e ){ c < m cos mX + S (m sin mX) (1) 

f (=o\ r J m= 0 

where r is the radial distance from the center of mass, 6 is the colatitude, A is the east 
longitude, and G - 6.673 x 10'" m 3 kg' 1 s' 2 is the universal constant of gravitation, while M E 
= 5.974 x 10 24 kg is the mass and R E =6.378 x 10 6 m is the equatorial radius of the Earth 
(9). The P f m (cos 6) are the associated Legendre polynomials, while the C tm and Sf,n are 
unnormalized spherical harmonic coefficients. Equation (1) can be more compactly written 

v = (2) 

r r J 

where Y (m] (d, A) = P lm ( cos 6 ) cos mX A) = R^ m (cos 6) sin mX , and 


'fm\ 


(2-<5 0m )(^-m)! Jp(r,0,A)/ T toll (0,A ) dv 

(£ + m)\ 


M E R e E 


( 3 ) 


with C lml = C tm and C fm2 = S fm . Here p (r,0,A) is the density and dv is an element of 
volume; the integral is over the whole Earth. The spherical harmonics are normalized 
according to the following equation: 


Rubincam 


3/27/01 


4 


jlY, m (6.M?dA 


47r(£ + m)l 


(4) 


where dA is an element of area on the unit sphere. Likewise, by analogy with (2) the 
gravitational potential from the core alone can be written 

_ Vjf y ((U) (5) 

r r ) 

where M c = 1 .947 x 10 24 kg is the mass and R c = 3.480 x 10 6 m is the radius of the core, 
with the K tmi being the spherical harmonic coefficients of the core. 

The density anomalies in the core will attract the density anomalies in the mantle, 
producing a torque L on the mantle. An element of torque will be given by 

dL = p(r,6,X)dv( rxVV c ) (6) 


where r is a point in the mantle and 

r = r(sin 0 cos A x + sin 9 sin A y + cos 6 i ) 

Here x, y, and z are unit vectors in a right-handed coordinate system fixed in the mantle, 
with z lying along the Earth’s rotation axis and x lying in the equator at zero longitude. 
For the first few harmonics the expression r x VV C can easily be found from the Cartesian 
formulation of the unnormalized associated Legendre polynomials. Writing F 2OI (0,A) = 
(2z 2 - jc 2 - y 2 )/ 2, etc., where r 2 = x 2 + y 2 + z 2 , results in such expressions as 
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etc., so that rxVV c can be written 


r x VV C 


GMc y f RA 

Rc r J 


i) 


(7) 


where the L) mi , U (mj , and L) m are coefficients given in Table 1. Integrating (6) gives the the 
total torque on the mantle: 


1 ■ r 2 sin 6 dXdddr (8) 


The density can also be expressed in terms of spherical harmonics: 


p(r,e,X) = ),A) 

l mi 


(9) 


By (7)-(9) the component of the torque along the rotation axis L z can be written 


L = 


4 7r{£ + m)\ 


(2-S 0 ,)(2(+l)(e-m)< 


Z^, G McRcl^ CS'Mf ' * 


tmi 


Rr 


l - I 


( 10 ) 


with similar expressions for L x and L v . 
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The nonsecular acceleration will be examined first. The time rate of change in the 
mantles’s rotation rate Q will be given by 


0,6 = L z ( 11 ) 

where C M is the moment of inertia of the mantle and currently Q = 7.29 x 10' 5 s' 1 . 

What is needed now is an estimate of the integral in (10) to use in (11). I will 
assume here that « |C, mi | for all coefficients (except perhaps the / = 1) terms, so that 

the C (mi coefficients refer to the mantle alone with a high degree of approximation. 
Moreover, assuming all the mantle density anomalies lie at the Earth’s surface probably 
tends to minimize the torque. Hence from (3) and (9) 

4 tc\ Rf p tmi (r) r ul dr 

C = — ^ (12 

(2 1 + 1 )M e R[ 


and setting r = R E in the above equation and (10) yield 


GM C M E -y (l + tn)\ z _ 


V R e j 


(13) 


as a conservative estimate of the torque by the core density anomalies on the mantle 
anomalies. The largest C (mi are the second degree terms. Setting /= 2 in (12), writing M c 

~ M e /3, and 5(R C /R E ) 2 ~ 1.5 in that equation along with C M ~ M^/3 in (1 1) give, with 
the aid of the entries in Table 1, 
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Q - 7.2— ^(K ln C 21l -K m C lu ) (14) 

as the approximate nonsecular acceleration of the Earth. The coefficients are known to be 
C 221 = 1.57 x 10' 6 and C 222 = -0.9 x 10' 6 (70); since only the order-of-magnitude is 
important here, they will be taken to be 10' 6 . The coefficients K 22l and K 222 are unknown as 
to sign, but the magnitude of the density anomalies in the core are expected to be on the 
order of Ap/p = 10' 8 -10' 9 , assuming the anomalies are spread evenly throughout the core 
(11). Using the smaller value in the equation analogous to (12) gives core coefficients K 2mi 

on the order of 10 10 . Plugging such values in (14) yield values of £2 = 7 x 10' 22 s' 2 , which 
is the observed order-of-magnitude for the nonsecular acceleration. If characteristic core 

flow is 10" 4 m s" 1 (7), and the blobs dissipate after traveling across the core (about 3000 
km), then the intrinsic anomaly timescale of the core variable field will be -1000 years. 

Suppose that the core gravitational field is responsible for the variable nontidal 
acceleration of the Earth’s rotation and that this acceleration acts on the observed 1000 year 
time scale. A typical change in mantle rotation speed will be the acceleration x time. The 
acceleration has been positive (2) for about the last 1000 years (see Figure 1). so that the 
change works out to be (ignoring other torques) about +0.03° per year in recent times. By 
action-reaction the torque will slow down the core, which will respond 10 times more 
quickly to the torque than the mantle, due to the core’s smaller moment of inertia (assuming 
it behaves like a rigid body), and in a sense opposite to the mantle. Hence one would expect 
the core to drift westward relative to the mantle at a typical speed of 10 times faster, which 
works out to be 0.3° per year. The observed drift of the magnetic field is westward at about 
0.2° per year (12), in good agreement with this mechanism. Moreover, the sign and rate of 
drift can be expected to vary; this is also observed (12) and occurs on the 1000 year time 
scale, in rough agreement with the nonsecular change in Earth rotation rate. (The speed at 
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which the core’s lumpy gravitational field washes over the mantle depends on the strength 
of the coupling and introduces another timescale apart from the intrinsic anomaly timescale 
over which the torque varies, but at 0.3° per year this is also happens to be about 1000 
years.) 

The core-mantle gravitational torque will also affect the obliquity e, which is the tilt 
of the Earth’s equator with respect to the ecliptic. Here the biggest contributor to the torque 
will come from the core field acting on the rotational flattening of the mantle, embodied by 
the J 2 =~ C 2 01 = 1 x 1 0' 3 coefficient in the external field, which is about a thousand times 
bigger than any other terms in the field ( 10). In this case the rate of change of the mantle’s 
obliquity with time will be given approximately by 


= LMLl 
CD. 


(^ 2 . 2 * - K my) 


which amounts to about 0.5 m yr' 1 . The present rate due to the Earth’s precession and 
orbital perturbations is about 14.7 m yr' 1 ( 13), so that core-mantle coupling accounts for 
only about 3% and does no great violence to the expected total. Being an internal torque, 
the rate of change of the obliquity of the core will be about ten times faster (about 5 m yr' 1 ) 
because of its smaller moment of inertia, and in the opposite direction. 

The core-mantle gravitational torque also produces polar wander. The magnitude of 
this wander is governed by the Liouville equations ( 14): 


m x + <J r m y 


<7 L 

r x 

n 2 ( c-a ) 


(15) 


a 2 ( c-a ) 


m y - <J r m . 


(16) 
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m x and m y are the x- and ^-components of the spin vector divided by Q, and where 


<* r 


C-A 


q. 


These equations ignore all effects on the pole except the torques. Here A is the Earth’s 
moment of inertia measured about an axis in the equator, and (C - A) I A = 1/300, so that a r is 
a frequency with a period of about 1 year. Differentiating (15) with respect to time and 
using (16) in the resulting equation gives 


m x + cr r 2 m x 


A. 

AQ 


+ <x„ 


AQ. 


(17) 


To make this equation easy to solve, I assume that the torque has the form L x - l? x e~ ai with 
a similar expression for L y , so that the torque varies on a characteristic timescale of a' 1 . 
Since this timescale is assumed to be hundreds of years or more, a«c r . Thus the first term 
on the right side of (17) can be ignored compared to the second. Assuming a particular 
solution m x = m x e~ a ' gives approximately 


m r 


<7 r AQ 


with a similar expression for the y-component. Differentiation the particular solution and 
using the expressions analogous to (13) yield a rate of polar wander of 
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■ A .A GC Z / 1 ' -j/ A r . A . 

m,x + m v y - — 7— f7 2 (^ 212 x - ^ 2ll y) 

C7 r i2/? £ 

which is once again dominated by the rotational flattening term J 2 . For K 1U ~ 10' 10 and 
(a/G f ) ~ 10' 3 the rate of polar wander is only about 0.05 cm yr ', which is well below the 
observed 10 cm yr' 1 thought to be induced by postglacial rebound. 

Gravitational core-mantle coupling appears to be a candidate for explaining the 
fluctuations observed on the 1000 year timescale in the Earth’s rotation rate. It may also 
explain the westward drift of the core-relative to the mantle, and would add a stochastic 
component to the rate of change of the mantle’s obliquity at perhaps the 3% level. 

However, this mechanism would not be able to explain any substantial north-south drift of 
the magnetic pole relative to the mantle, since this is given by the change in the core’s 
obliquity relative to the mantle, which would only be about 5 m yr 1 . The contribution to the 
presently observed polar wander is tiny, only about 0.5%. 

The above calculations assume that the viscous drag of the outer core on the mantle 
is negligible. The actual viscosity is unknown. The characteristic spin-up time for a core 
rotating at a speed different from the mantle is between 10 years and 10 5 years for a 
plausible range of viscosities (15). Thus only for the lower estimates of viscosity is its 
neglect justified. 

The calculations also assume that the outer core rotates like a rigid body relative to 
the mantle, when in fact it is a convecting liquid. Further, the convecting blobs in the outer 
core will exert pressure and electromagnetic as well as viscous forces at the core-mantle 
boundary (5,16)\ these are not including here. The sum total of these effects represents 
future avenues of research, as do the consequences of this coupling for the long-term 
obliquity of the Earth (4,17) and the rotation of the solid inner core (18, 19). 
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Table 1. The coupling coefficients L x (mi , U tmi , U fmi of equation (7) for 1 = 2. 
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“ 2 Ar 22 2 
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Figure 1 

Nonsecular change in length-of-day (ALOD) in milliseconds vs. years AD, derived 
from subtracting a secular trend of +1.6 ms cy' 1 due to a combination of tidal friction and 
(possibly) postglacial rebound from the total observed ALOD. Based on the work by F. R. 
Stephenson and L. V. Morrison ( 2 ). 

Figure 2 

Schematic of the density anomalies in the liquid outer core and the anomalies in the 
solid mantle, looking down from the north pole. The observed positive acceleration of the 
mantle over the last millenium (outer counterclockwise arrow) and the current westward drift 
of the core (inner clockwise arrow) may be due to internal gravitational torques between the 


mantle and core. 
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Gravitational Core-Mantle Coupling and the Acceleration of the Earth 

POPULAR SUMMARY 

The tides are slowing down the rotation of the Earth, overwhelming a smaller 
positive acceleration from postglacial rebound. According to historical records such as 
ancient and medieval observations of solar eclipses, the signal from these two secular effects 
appears to be significantly corrupted by a source (or sources) of an episodic or perhaps 
periodic nature operating on a 1000 year timescale. This nonsecular source causes 
accelerations which are the same order-of-magnitude as the tides, ~6 xlO' 22 s' 2 . Coupling 
between the liquid outer core of the Earth and the mantle has long been a suspected reason 
for changes in the Iength-of-day. Research heretofore has been aimed at viscous, 
electromagnetic, and topographic coupling at the core-mantle boundary, mostly on the 
decadal timescale, and at gravitational coupling not related to the problem treated here. 

The present investigation focusses on the gravitational coupling between the density 
anomalies in the convecting liquid outer core and those in the mantle and crust as a possible 
cause for the observed nonsecular acceleration on the millenial timescale. (Hereafter I use 
"mantle" as a shorthand for "mantle and crust.") The basic idea is as follows. There are 
density inhomogeneities caused by blobs circulating in the outer core like the blobs in a lava 
lamp; thus the outer core's gravitational field is not featureless. Moreover, these blobs will 
form and dissipate somewhat randomly. Thus there will be a time variability to the field. 
These density inhomogeneities will gravitationally attract the density anomalies in the 
mantle. 

The outer core and mantle may not be in a state of gravitational equilibrium if the 
blobs form and dissipate faster than either the outer core or mantle can react to achieve a 
state of zero torque. In other words, the mantle and outer core may always be in a state of 
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disequilibrium due to the stochastic nature of "core weather." The outer core and mantle 
will try to align with each other to make the torque between them zero; but before this can 
happen, the gravitational field of the outer core may change, so that the core and mantle 
never catch up with each other. Hence it may be that the outer core is forever torquing the 
mantle, changing its rotation state and perhaps producing the observed nonsecular 
acceleration and changing the length-of-day. By action-reaction the outer core will also 
change its rotation state. 

A quantitative calculation indicates that the gravitational torque on the mantle from 
the liquid core may be sufficient to explain the accelerations and decelerations in the Earth’s 
rotation on the 1000 year timescale. The corresponding torque by the mantle on the core 
may also explain the westward drift of the magnetic field of 0.2° per year. Gravitational 
core-mantle coupling would stochastically affect the rate of change of the Earth’s tilt by just 
a few per cent. Its contribution to polar wander would only be about 0.5% the presently 
observed rate of 10 centimeters per year. 


